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1. Introduction 

The strong sector of the Standard Model Lagrangian contains a number of inputs that are 
a priori unknown and must be determined from experiment. For these fundamental constants; 
the quark masses and the strong coupling constant, our knowledge is currently rather imprecise. 
Their determination is complicated by confinement in QCD, so the quarks and gluons cannot be 
observed as isolated particles. We can only determine these fundamental parameters by solving 
QCD for observable quantities such as hadron masses, as functions of the quark masses and the 
coupling (or alternatively, the lattice spacing). The High Precision QCD (HPQCD) collaboration 
has a program for calculating the values of these parameters using the numerical techniques of 
lattice QCD simulations combined with multi-loop perturbative renormalisation. Previously the 
masses have been determined at the one-loop level |[T[ ^ ||], and the strong coupling constant was 
determined from one observable in a two-loop calculation [Q]. Recently the determination of the 
strong coupling was improved with three-loop perturbative matching and with input from 28 lattice 
observables in simulations at three different lattice spacings, resulting in an accuracy of just over 
1% lH |]. This writeup covers progress on the determination of the light quai^k masses, where we 
push the perturbative matching calculation to two loops, that is, next-to-next-to-leading order. 

Precise knowledge of quark masses constrains Beyond the Standard Model scenarios as well 
as providing input for phenomenological calculations of Standard Model physics. The strange 
quark mass, in particular, is needed for various phenomenological studies, including the important 
CP- violating quantity e'/e [0], where its uncertainty severely limits the theoretical precision. 

Previously, shortcomings in the formulation of QCD on the lattice and limitations in comput- 
ing power have meant that lattice calculations were forced to work with an unrealistic QCD vac- 
uum that either ignored dynamical (sea) quarks or included only u and d quarks with masses much 
heavier than in nature. This condemned determinations of most phenomenologically-important 
quantities, including the quark masses, to rather large systematic errors (10-20%) arising from the 
inconsistency of comparing such um^ealistic theories with the necessary experimental input. The 
determination presented here uses simulations with the improved staggered quark formalism that 
generates a much more realistic QCD vacuum with two light dynamical quarks and one strange 
dynamical quark. Staggered quarks are fast to simulate. They keep a remnant of chiral symmetry 
on the lattice, which prevents the occurrence of exceptional configurations, and which allows sim- 
ulations at much smaller quark masses. The bare quark masses in the simulations were fixed using 
chiral perturbation theory to reliably extrapolate lattice results to the physical point [[T]]. Working 
in the region of dynamical u/d quark masses below ms/2 and down to m.v/8 gives control of chiral 
extrapolations and avoids the large systematic errors from dynamical quark mass and unquenching 
effects that afflicted previous studies using other lattice discretisations. 

The dominant systematic eiTor in the determination of the MS masses in [|l|] came from un- 
known second- and higher-orders in the perturbative matching. Some progress was reported on the 
chiral fits at the lattice meeting ||8|], however that analysis is not used here; we continue to employ 
the bare quark masses given in [|T]]. Significant progress on the reduction of the systematic eiTors 
is reported here, due to our computation of the second-order perturbative matching coefficient, the 
first determination at this order of a "kinetic" quark mass in any lattice theory (the zero-point ad- 
ditive renormalisation for Wilson fermions was previously determined at two-loops in [^, and for 
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Static quarks in [p^). 

The staggered quark formalism does present several challenges, which have been tamed with 
an aggressive program of perturbative improvement. With the naive staggered action large dis- 
cretisation eiTors appear, although they are formally only 0'{a^) or higher {a is the lattice spacing). 
In the case of the unimproved staggered action the renormalisation of lattice operators to match 
to continuum quantities were also frequently large and poorly convergent in perturbation theory. 
This was true, for example, for the mass renormalisation that is needed here. It turns out that both 
problems have the same source, a particular form of discretisation error in the action, called "taste 
violation," and both are ameliorated by use of the improved staggered formalism [11]. The per- 



turbation theory then shows small renormalisations [ |12| , |13| , |14J, [15|] and discretisation errors are 
much reduced [^, [l7|, 18 1. Empirically, taste violation remains the most important discretisation 
error in the improved theory, despite being subleading to "generic" discretisation errors. The Gold- 
stone meson masses we will discuss here are affected by this at one-loop in the chiral expansion. 
Staggered chiral perturbation theory (SZPT) [jl9|, ^ ^] allows us to control these effects and 
reduce discretisation errors significantly. 

A potentially more fundamental concern about staggered fermions relates to the need to take 
the fourth root of the quark determinant, in order to convert the four-fold duplication of "tastes" 
into one quark flavor. One might imagine that the fourth root introduces nonlocalities which pre- 
vent decoupling of the ultraviolet modes of the theory in the continuum limit. However, evidence 
is amassing that demonstrates that the properties of the staggered theory, with the fourth root, 
are equivalent to a one-flavor theory, up to the expected discretisation eiTors. These are due to 
short-distance taste-changing interactions, which are mediated by high-momentum gluons [11] 
(the locality of the free-field staggered theory is trivial, and is made manifest in the "naive" basis 
used in [|n|]). One should not be surprised that nonlocalities do not arise, precisely because the 
staggered quark matrix is diagonal in the taste basis, up to those small, short-distance (and cal- 
culable) corrections. It has been demonstrated that perturbative improvement of staggered actions 
correlates exceedingly well with non-perturbatively measured properties of the staggered fermion 
matrix, providing clear support for the correctness of the fourth-root procedure. This includes the 



measured pattern of low-lying eigenvalues of the staggered matrix [ ]23| , |24[ , |25| ], and the measured 
pattern of taste-violating mass differences in the non-chiral pions [^]. 

The rest of this paper is organised as follows. The methodology of the calculation is discussed 
in the following section, including a general discussion of how to obtain the matching factor which 
connects the bare quark mass mo(a) to the MS mass ni^^{}x), using the pole mass as an intermediate 
quantity. We also derive the two-loop anomalous dimension for the bare quark mass in the lattice- 
regularised theory, and give complete results for the two-loop matching coefficients. Section || 
gives preliminary results for the light quark masses, with a preliminary analysis of the systematic 
uncertainties, including a technique to make a rough estimate of the third-order perturbative cor- 
rection. Section ^ compares our results with other recent determinations of the strange-quark mass. 
An appendix provides some additional information concerning the evaluation of the multi-loop 
diagrams, including explicit expressions for the two-loop quark mass renormalisation in terms of 
the IPI self-energy, a discussion of the techniques used for generating and evaluating the neces- 
sary multi-loop integrands, and some detailed numerical results which provide an indication of the 
many consistency checks that we have appUed to our results. 
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2. Structure of Calculation 



The quark masses ai^e not physically measurable, and as such are only well-defined in certain 
renormalisation schemes, such as the MS mass m^^{}i.), evaluated at some convenient scale jX. The 
light quark MS masses are determined here by multiplying the chirally-extrapolated bare masses 
in lattice units, amo, by the inverse lattice spacing a, and by the appropriate perturbative matching 
factor Z„j{iJ.a,moa), which we compute to two-loop order: 

{mo a) 



m 



MS 



(M) 



Zm{}la,mQa), 



(2.1) 



where the bare mass mo (a) is cutoff dependent. The bare masses used here were determined in 
an extensive chiral perturbation theory analysis of the MILC Asqtad data that was discussed in [|l|. 



27| , [28| , |29[ ]. These bai^e masses were previously used with the one-loop perturbation theory result 
for Zm to extract the MS masses in [|l|] ; the final result for the strange quark mass reported there 
was 76(0)(3)(7)(0) MeV where the respective errors are from: statistics; simulations systematics 
of which the most important are chiral fitting and lattice spacing; an estimate of the unknown 
two-loop perturbative errors and an estimate of the uncertainty due to electromagnetic and isospin 
contributions to the pion and kaon. That the eiTor coming from the lattice spacing a is so small 
is a distinguishing feature of this calculation. The lattice spacing a is one of the five simulation 
parameters, and an important one because it sets the simulation's mass scale. In our eai^lier light 
quark masses and a, analyses, we set the lattice spacing by comparing a simulated T mass splitting 
{e.g., T' — T) with experiment. Here we continue this practice, although the lattice spacings derived 
from our T splitting have been shown to agree with spacings derived from a wide variety of other 
physical quantities: ten in all, including the pion and kaon leptonic decay constants, the Bs mass, 
and the ^l baryon mass iQ |3^]. All of these different quantities agree on the lattice spacing to 
within 1.5-3%. 

The largest error in our previous determination of the quark masses [|l|] was from the pertur- 
bative matching. Here that is addressed by the calculation of Z,„ at two loops. We do this in two 
stages, using the pole mass M as a matching quantity to connect the lattice- and MS-regularisation 
schemes. We also use our previous determination of the relation between the lattice bare coupling 
and the renormahsed couphng av{q*), defined by the static potential, to reorganise both sides of 
the matching equation into series in terms of (Xv{q*) at an appropriately determined scale. 

We begin by recalling the relation between the MS mass and the pole mass M, which is known 



through three loops 



|, |34], |35|]. We require it to second order, a result that was first 



obtained in |32] (expressions for the relation at arbitrary /i are conveniently given in [34]) 



71 



+ . 



(2.2) 



where the one- and two-loop coefficient functions z\{lx/M) and ziijx/M) are reduced to a set 
of terms with different colour structures [in the following Cf = (A'^^ — \ )/{lNc), Ca = Nc, and 

r = 1 /2] 



zi = CfZf 

Z2 = CpZFF + CfCaZfa + CFTneZFL + CfTzfh, 



(2.3) 
(2.4) 
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and where the contribution zfh from an internal quark loop with the same flavor as the valence 
quark is spht off from the contribution zfl of internal quark loops with different flavor (these are 
taken here to be degenerate in mass, though this is easily generalised). The total number of flavors 
is «y = + 1. The individual functions are given by 



Zf 
zff 
Zfa 



1 

128 



1111 I 1 r I 3 r 3 



384 ' ~'~ 2 ^2 + g ^3 ■ 



:C2log2- 



185/ „_il/2 
95 <-nM 



Zfh 



vi + i^2-2A(l) + i|V + ^^' 



(2.5) 
(2.6) 
(2.7) 
(2.8) 
(2.9) 



where i^M = log(At^/M^), and where the function A(r) gives the dependence of the renormalisation 
factors Zfl and zfh on the quark mass in an internal fermion loop (sea and valence, respectively), 
with Tsea = wisea/"^vaience ■ An cxact integral expression for A(r) can be found in [32] and it is plotted 
in the appendix. 

Next we outline the computation of the relation between the pole mass M and the bare quark 
mass mo (a) in the lattice-regularised theory, which has the form 



M = mo [ 1 + ttL (A 11 log mod + A lo) + (A22 log^ mo^ + A2 1 log moa +^20)+- 



(2.10) 



The coefficients of the logarithms are determined by the two-loop anomalous dimension in mo (a), 
which in turn can be determined from the known anomalous dimension of the MS mass, as de- 
scribed below. We have previously computed the one-loop matching factor, with the result ^ [I^] 
Aio ~ 0.5432 (neglecting corrections of C?((a'mo)^)). What is new in this paper is the determination 
of the two-loop term ^20- 



As with the two-loop continuum matching factor in (2.2), A20 depends on the number of quark 
flavors in the sea, and on the ratios of the sea quark masses mgga to the valence quark mass myaience- 
However, as we demonstrate explicitly in the appendix, the mass dependence in A20 cancels pre- 



cisely in the matching to the continuum relation ( |2.2[ ) for mjea/mvaience in the range necessary. 
This follows from the fact that, in the limit where the energy scales are large, the net renormal- 
isation factor connecting mo(fl') to m'^^(/i) probes internal loops at scales large compared to the 
internal quark masses; hence the dependence on quark masses of the intermediate renormalisa- 
tion factors connecting to the pole mass M are infrared effects that are identical in the lattice- and 
MS -regularised theories. 



Before we give our results for A20, and the final matching factor Z,„ in (2. 1 ), we consider how 



to determine the coefficients of the logarithmic terms in ( |2. lOj ). The pole mass is an RD invariant 
(equivalently does not depend on a) so taking logs and differentiating with respect to log a, and 
neglecting any scale dependence in the coefficients A„,„ (which would arise as discretisation errors 
that go as powers of moa < 0.05, which are neghgible for our purposes), one has 



log mo 

dloga dloga 



+ ■ 



(2.11) 
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where the square bracket is from (2.1C), and the first term above can be identified with the anoma- 
lous dimension equation: 

^ 70 aL{a) + y',al{a) + ^(a|) (2.12) 

The L's indicate the lattice-regularisation scheme. The sign comes from differentiating with re- 
spect to the logarithm of the lattice scale a^^ (momentum units), instead of pL. Only the leading 
term in the anomalous dimension is scheme independent, 70 = 3Cf/{2n). Comparison of the two 
preceding equations leads to the identification 

All = -70, A22 = ^Ali-An^, A21 = -7^ +A?i -Aio -AnY (2.13) 



where the jS -function arises from an implicit derivative of the bare lattice coupling in ( 2.1 1| ) (j3( 



11 — !«/). It remains to determine "j^. This can be found from the known MS anomalous dimen- 
sion. 



making the substitution = a \ and using the relations mj^(a; ^) =mo{\+CmOc), and a-^{a ^] 



aL{ \ +CaOCL). By comparison to (2.12) we obtain 



rf = rr'+c„7o-c,„|^. (2.15) 

271 

(2) 

The one-loop renormalisation constants C„, =Zm (jLt = l/a) =0.1188 andCa =4.753 — 0.3316?!/^ 
for the improved gluon and staggered quark actions have been previously calculated [jl], [T2|,|5[,[6|]. 
Knowing the logarithmic terms due to the lattice anomalous dimension provides a useful cross- 



check on our evaluation of the two-loop renormalisation factor (2.1C). We compute the two-loop 
lattice diagrams as a function of moa, and subtract the known logarithms, in order to isolate the 
remaining term A20, which must be finite as niQa — > 0. Additional checks are provided for diagrams 
which have a leading log^(a;nio) term, which arises from the infrared limit of both the outer and 
inner loop integrals, and which is therefore an infrared quantity; the coefficients of the double loga- 
rithms in the individual diagrams are available in Feynman gauge from the original MS calculation 
of Tarrach and we have verified that these are reproduced in our lattice calculation. 

The necessary two-loop diagrams are shown in figure [l[ Further details on our evaluation of 
these diagram are given in the appendix. We can easily evaluate them for various gluon and quark 
actions, using our automated methods for generating the lattice Feynman rules. In this paper we 
give results for the Symanzik improved gluon with improved staggered quarks (Asqtad) and SU(3). 
Our result for the matching term A20 is, in the limit of vanishing sea quark mass, 

A20 = 6.09 -0.15?if- 0.03, (2.16) 

where the last term corresponds to an internal quark loop containing one (massive) valence quark 
flavor. The uncertainties arise from a numerical evaluation of the two-loop integrals. 
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D6 



e 



D7-11 D12 D13 



A XX 



D2: 



D14-18 D19 D20 

D21 D22 D24 

j D26 D27 *l j» D28 



CTl CT2 CT3 CT4 CT5 CT6 

Figure 1: The two-loop diagrams that contribute to the two-loop mass renormalisation. The numbering is 
consistent with [^. Blue lines are valence quarks (on-shell externally), green are gluons and red are sea 
quarks. The large brown vertices are stand-ins for the five one-loop gluon propagator diagrams (gluon and 
ghost bubbles and tadpoles and the measure term). We evaluate all five simultaneously with an internal 
subtraction. The crosses represent interaction vertices generated by perturbative expansion of tadpole and 
other renormalisation factors in the gluon and quark actions. 



When the lattice renormalisation factor (2.10) connecting the bare mass to the pole mass is 



combined with the equivalent continuum expression ( [2.2| ) for the connection to the MS mass, the 
logs of mo drop out, as expected, since these are infrared effects that are identical in the intermediate 
lattice and continuum matchings. We also reorganise the couplings to the ay scheme at some scale 
q*, where q*a is a function of iJ.a that is determined according to the BLM scheme [|39|]. This 
leaves an expression with logarithms only of jxa and q*a. The final expression for the perturbative 
matching factor is then: 

Z,n{lJ.a,moa) = 1 +zi?^(/la)ay (<7*(jUa)) +Z,lf'(/la)ay + ^((moa)'^, (moa;)^a^, a^), (2.17) 



011/7 



High Precision Fundamental Constants using Lattice Perturbation Theory 



Quentin Mason 



where Z,n'{na) = 0.1188(1) — 2/nlog{na) + ^{moa) was derived previously, [p, and the new 

(4) 

information presented here is the expression for Z,„ 

Zif'(jua) =Z2log2(jUa)+Zilog(Ata)+Zo, (2.18) 

where the latter coefficients depend on logarithms of q*a 

Z2 = 0.76-0.034«/ (2.19) 
Zi = -0.40- 0.028n/ + (- 1.12 + 0.0675?!/) lng*a (2.20) 
Zo = 2.09 - 0.014?!/ + (0.21 - 0.013?i/)ln^*a. (2.21) 

We have computed the optimal choice for q*a in the BLM scheme, as a function of jXa, using an 
exact evaluation of the average momentum scales in the continuum self-energy diagrams, according 
to [ ^ ] (this is an improvement over the inexact calculation in [|]]), and a numerical evaluation of 
the average scales in the lattice self-energy (unchanged). A typical value for the matching scale is 
q*a = 1.88 at jxa = 1. 



3. Results 



The bare lattice masses for the strange and up/down quarks, on the MILC "coarse" and "fine" 
lattices, are given in For the strange quark, these are amos = 0.0390(1) (20) /mqo and amos = 
0.0272(1) (12)/mo/, on the coarse and fine lattices respectively, where mqc = 0.85488 and mq/ = 
0.86774 are tadpole normalisation factors. The uncertainties are lattice statistical and systematic 
errors, respectively, the latter due mainly to chiral extrapolation/interpolation. The lattice spacings 
can be found in = 1.596(30) GeV, and a^^^ = 2.258(32) GeV. 

Following conventional practice, we quote the light quark MS masses at the scale n =2 GeV, 



taking three active flavors of quarks (?i/ = 3). This results in two-loop coefficients in ( |2.17[ ) 



of Zm I coarse — 

1.939(4), and Zm^lfine — 2.270(4). Finally, we require the couplings at the rel- 
evant scales, and for this purpose we use the recently determined value ay"'^~^^(7.5 GeV) = 
0.2082(40) [gl. We find av{c]tom-sc) = 0.2925(92) and aviql^,) = 0.2713(76). 

Putting all of this together, we obtain the following preliminary values for the MS strange- 
quark mass, using the coarse and fine lattices as input 

?wf^(2GeV) = 83(5) MeV [coarse], ?72p(2 GeV) = 85(4) MeV [fine], (3.1) 

where the errors here are just from the simulation systematics. 

Following [|T]], we consider continuum extrapolations of the central values in (3.1) based on 



the form of the expected leading discretisation en^ors, which are of O(aya^), and compare with a 
extrapolation in ayO^. These two forms yield almost identical extrapolations, to a central value of 
87 MeV, with an error that is folded into the quoted lattice systematic errors above. We are currently 
considering refinements to our estimates of the continuum extrapolation and the incorporation of 
third order perturbative eiTors. For this preliminary result we quote as our continuum extrapolation 

?7ip(2 GeV) = 87(0) (4) (4) (0) MeV (3.2) 
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where, following [[T]], the respective eiTors are statistical, lattice systematic, perturbative, and elec- 
tromagnetic. We have assigned a relative eiTor from the estimated third-order perturbative matching 
of roughly ±2ay{q*) ^ 5%, down from the 10% perturbation theory eiTor in the one-loop determi- 
nation in [|l|]. We anticipate to reduce somewhat the perturbation theory eiTor with a more careful 
study of the systematics of the jx dependence. 

Our result for the ratio of strange to up/down quark masses is unchanged at 27.4(1)(4)(0)(1) 
from [I], since the renormalisation factor is mass-independent, as we have verified above explicitly 



through two-loops. Using the value mu/m^/ = 0.43(0)(8)(1)(0) from [28| this gives: 

mf^(2 GeV) = 1.9(0)(1)(1)(2) MeV 
m^(2 GeV) =4.4(0) (2) (2) (2) MeV 



4. Discussion and Conclusions 



Perturbation theory has once again shown itself to be an essential tool in high precision phe- 
nomenological calculations from the lattice. The two-loop lattice diagrams were conquered with 
a combination of algebraic and numerical techniques in this first-ever determination of the multi- 
plicative two-loop pole mass on the lattice. When combined with the known continuum matching 
from the pole mass to the MS mass, a very accurate determination of the light quark masses was 
possible. The results presented here have a number of distinguishing features: two-loop perturba- 
tion theory, ?iy = 2 + 1 simulations with two degenerate light quarks and a heavier strange quark, 
very small light quark masses from m^ /S to nis/l which enabled a partially quenched chiral fit with 
many terms to thousands of configurations, and extremely accurate determinations of the lattice 
spacings, which are equal within the errors when set from any of 10 different quantities, ranging 
from the very lightest states all the way up to heavy mesons and baryons [Q, 

Most notable amongst our results is our new value for strange quark mass m^^(2 GeV) = 
87(0)(4)(4)(0) MeV, where the respective eiTors are lattice statistical, lattice systematic (mostly 
due to the chiral extrapolation/interpolation), perturbative, and due to electromagnetic effects. The 
two-loop matching has increased the central value with respect to the previous determination in [|l]] 
by about two standard deviations, based on the previous estimate of the perturbation theory un- 
certainty (which was ±7 MeV). We believe the present estimate of the perturbative uncertainty of 
±4 MeV can be reduced somewhat by estimating the third-order perturbative correction, along the 
lines that we have already explored, in a preliminary way, here. 

The strange quark mass determination has historically generated some controversy, with some- 
what different values being obtained from different approaches. An obvious advantage of our result 
is that it has been obtained with the correct description of the sea, that is, with ?iy = 2 + 1 flavors 
of dynamical quarks. There is only one other result with the correct number of flavors in the sea, 
which is due to the CP-PACS and JLQCD collaborations, who reported a value at this conference 
of m^^(2 GeV) = 87(6) MeV [p0|], although the enw was not very well quantified, and does not 
appear to include an estimate of corrections due to two-loop and higher-orders in the perturbative 
matching. 

It appeai^s the the strange-quark mass extracted from simulations with only Uf = 2 flavors of 
sea quarks are systematically higher than the estimates with the correct rif (noting that the two- 
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flavor determinations were also done with nonperturbative definitions of the quark mass), although 
the other systematic enws are too large to allow for a definitive assessment. The two-flavor de- 
termination from the QCDSF-UKQCD collaboration, is mp(2 GeV) = 100-130 MeV Q, the 
ALPHA collaboration value is 97(22) MeV [||], and the Rome value is 101(8)(tf ) MeV [||]. 

The next major evolution in our lattice determination of the light quark masses will come 
from using a third lattice spacing, the "super-fine" configurations already partially implemented by 
MILC, and planned by UKQCD. This will hopefully reduce the size of the chiral and systematic 
errors due to taste-breaking, and vastly improve the quality of the continuum extrapolation. The 
HPQCD collaboration plans to continue its accurate determinations of the fundamental constants 
by generalising this calculation in the first instance to heavy quark masses to complete the strong 
sector of the Standard Model Lagrangian with two-loop calculations of all the quark masses with 
the already finished three-loop strong coupling constant. Further generalisation by insertion of 
appropriate operators will give the two-loop vector and axial-vector renormalisations, necessary for 
improving the accuracy of the decay constants fs, fo etc and associated form-factors. Eventually 
HPQCD plans therefore to have two-loop accurate determinations of many of the CKM matrix 
elements in the next few yeai^s. 
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A. Perturbative Structure 



We present some details of the lattice size of the pole mass renormalisation factor, ( |2. 10[ ). We 
define the full lattice quark propagator to be: 

G{p,moy^ =if+mQ + 'Ltot{p), (A.l) 

where Ztot is minus the usual IPI truncated two-point function. The lattice dispersion relation im- 
pUed by p will be kept as a general function of the lattice momentum p; for unimproved staggered 
quarks, for instance, p^ = sm(p^). Make the spinor decomposition: 

Itotb) =Si(p) + 0?+mo)l2(/'), (A.2) 

where both £i and £2 are implicitly functions of mo, and both are Dirac scalars. Note that £2 
is only part of the wavefunction renormalisation. The pole mass is defined by the all-orders on- 
shell condition, corresponding to the location of the pole in the propagator. At tree-level this is 
"—if = wio" (a very common but abusive notation). As is conventional, we work at zero external 
three momentum, p = {ptS^), and use the positive energy spinor projector (Euclidean), (1 + 7f)/4. 



Then it is convenient to rearrange (A.l) and (A.2) to get 



nP.)^-ip.=mo + ^^^ (A.3) 
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This relation can also be applied to theories with an additive mass renormalisation, such as Wilson 
quarks, by absorbing the momentum dependence of the additive mass into P{pt)- 

We recursively solve for the on-shell pole mass M, defined by the renormalised energy 
at zero three-momentum 



pt = iM, where M = M^°^ + sIm'^^^ + g^M 



(4) 



(A.4) 



is the perturbative expansion, with a corresponding series for other quantities (including e.g. p) 
iM^^^ andP(/M(°^) =mo). 

We first consider the expansion of (|A.3|) through one-loop, which reads 



(0) 



Hence 



M 



PiiM^'HiglM('^)=mo + glzf\p^'^). 

dP{ix) 



(2) _ 7(0)^(2). (0) 



(0) 



7:^>l},\p^^^), whereZ^^ 



-1 



dx 



f=M(0) 



(A.5) 



(A.6) 



is the tree-level wave function residue (e.g. — 1/ cosh(M(°') for unimproved staggered quarks). 
To find the location of the pole at two-loops requires solving (AJ_) self-consistently 

P(/m(«) + /g2M(2) + iglM^') ) = ,„o + + glm(') , 

where the right-hand side above represents a self-consistent expansion of the right-hand side of ( |A.3| ). 
Part of the 0(go) ^^^^ arises from the one-loop piece of £1 when it is evaluated at the one-loop- 
corrected on-shell energy pt, determined by ( [/^.4[ ) and ( |A.6D . By Taylor expansion: 



zf{p^)=Yf{iMi^^)+iglM^'^^-^^ 



+ 



The result for the two-loop contribution to the pole mass is therefore 



.(0) 



;„(4) _ 1(^(2) )2^^^ 



(A.7) 



(A.8) 



(the second-term above is a correction of O((amo)^) for Asqtad), and 



m 



(4) 



zS')(/M(°))+£P(/M«")zi,°)<'Tr 



where the identity 



z(°)Tr 



l + y,a£|g(A 
4 dpt 



was used. In the continuum, Z^^ = 1 and ( |A.10| ) is the complete one-loop expression for the wave 
function renormalisation there. On the lattice, which has a tree level wavefunction renormalisation 
the result has some additional terms: 



\ + ytdZ^^{pt 
4 dpt 



OPt 



(A.9) 



p,=iMW 



(A. 10) 



7-1 



7(0)- 



l + y,^£p)(p) 
4 dpt 



+ 



(2) d^Pjix) 

dx^ 



P=P 



.(0) 
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Figure 2: Schematic representation of an IR subtraction, where appropriate traces of the self-energy with an 
energy-projector are impHcit in this schematic representation. 



Unfortunately ( |A.10| ) is IR divergent for both lattice and continuum. This infrared divergence 
precisely cancels against an IR divergence in the two-loop nested-rainbow diagrams (D21, D22 
and CT4 in figure This parallels the cancellation of an infrared divergence in the continuum 



two-loop diagram, D21, § § » > > which is likewise rendered finite by the iteration of 
the the one-loop self-energy, which generates a "counterterm" given by the one-loop mass shift 



times the one-loop wave function residue, analogous to the terms in (A.9) 



In this connection, we note that the continuum pole mass was shown to be infrared finite at 



two-loops by Tarrach 1 3 1 ] ; an all-orders proof of the finiteness of the on-shell self-energy has only 



recently been established, by Kronfeld [36]. 



On the lattice, the infrared cancellation in ( |A.9| ) has a schematic representation in diagram- 
matic form, as shown in figure ^. We numerically evaluate the integral for the two-loop diagram 
on the left in figure ^ with a subtraction, in its integrand, of the product of the independent inte- 
grands for the one-loop self-energy, and its derivative. This grouping, as indicated in figure is 
IR finite, and does not require any infrared regulator. We obtain a powerful check of this result 
by noting that this combination generates a leading logarithmic contribution to the anomalous di- 
mension of the mass which goes like log^(a^mQ), and whose coefficient is identical to that of the 
infrared-subtracted combination in the continuum, which can be found in Feynman gauge from the 



MS results in [31| 



All the diagrams for the two-loop self-energy, figure |T], were generated and evaluated inde- 
pendently by two of us. The Feynman rules for the highly-improved actions are exceedingly com- 



plicated, and were generated automatically using computer algebra-based codes [|14j]. Moreover, 
one of us has produced an algorithm which automatically generates the Feynman diagrams them- 
selves [^7|]. The two-loop integrals were evaluated numerically, using the adaptive Monte-Carlo 
sampling of method of VEGAS [^]. Practically, when performing the subtractions of figure ^ the 
integrals are much more convergent with a (4D) spherical transform and a further ^new = log|fc| 
transformation. The effect of the associated Jacobian's is to regulate the IR with k"' and a spherical 
cut-off at small \k\. We found that e^^" and smaller were sufficient to show the integrals inde- 
pendent of this cut-off. Almost all diagrams benefited from using the spherical transform, though 
some of the other "continuum" like diagrams are better behaved with the "log-spherical" transform 
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Asqtad multiplicative mass renormalisation at two-loops 



6.5 



6.4 



6.3 



6.2 



6.1 



5.9 



0.001 



fit to Ajo 
QM total 
HT total 



0.01 

Quark Mass [inQCi) 



0.1 



Figure 3: The total of the gluonic part of the two-loop contribution to the pole mass on the lattice with 
the Asqtad action, evaluated independently by two of the authors, in units of a^, for varying quark mass 
(mofl). The data is shown after subtracting the known logarithms: A22log^OTofl + A21 logmofl, and should 
therefore be the constant, A20, with no mass dependence. The absence of observed lattice artifacts of the 
form (m{)fl)"log'mofl, / < 2 justifies our assumption in the derivation that for small masses the lattice is a 
mass-independent renormalisation scheme. 



for very large numbers of integrand evaluations. We have checked that varying the cut-off makes 
no difference to the result. A powerful additional cross-check was provided by an explicit verifi- 
cation that our results are gauge independent, which we established numerically for two different 
bare-quark masses in three covariant gauges: Feynman, Landau and Yennie. 



We show results for the two-loop part of the pole mass, ( |2. lOD , coming from diagrams without 
fermion loops, in figure |[ As described in Sect. 2, we can test our calculation by subtracting the 
known logarithms in rtiQa, in order to expose the remaining factor A20, which must be finite in the 
limit niQa —>■ 0. Figure |3| shows results for the gluonic part of A20 over a wide range of bare masses, 
which clearly shows the expected limiting behaviour. 

A further stringent check of our evaluation of the diagrams with internal fermion loops (dia- 
grams 12, 13, 19 and 20 in figure |T]) is achieved by computing over a range of sea quark masses. As 
described in Sect. 3, the mass dependence in the intermediate renormalisation from the bare mass 
to the pole mass M should cancel against the renormalisation from M to the MS mass. We define 



A20(''sea) = A2o(0) + ^ Aiattice (''sea) , 



"^valence 



(A. 11) 



and compare with the analogous continuum function A(rsea)- We plot our results in figure ^, over a 
very wide range in r^^a, which shows the expected cancellation. 
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10 




10~' 10° 10' 10^ 



r 

Figure 4: Comparison of the sea loop-quark mass dependence of pole mass renormalisation matching fac- 
tors, on the lattice-side of the matching, and on the continuum side. The difference between the contin- 
uum squares and lattice crosses (calculated at a/«vaience — 0.001), is smaller than the errors. The variable 
= wJsea/wJvaience- The solid lines show limiting forms of the dilogarithmic but analytic continuum function. 
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